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Abstract 

We study the energy fructuations of a spatially homogeneous SU(2) Yang- 
Mills-Higgs system. In particular, we analyze the nearest-neighbour spacing 
distribution which shows a Wigner-Poisson transition by increasing the value 
of the Higgs field in the vacuum. This transition is a clear quantum signature 
of the classical chaos-order transition of the system. 
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1 Introduction 



In the last years there has been much interest in classical chaos in field 
theories. It is now well known that the spatially uniform limits of scalar 
electró dynamics and Yang-Mills theory exhibit classical chaotic motion 1 ^ 8 ). 
On the other hand, in field theories, less attention has been paid to quantum 
chaos, i.e. the study of properties of quantum systems which are classically 
chaotic 9 -*. 

The energy fluctuation properties of systems with underlying classical 
chaotic behaviour and time-reversal symmetry agree with the predictions of 
the Gaussian Orthogonal Ensemble (GOE) of random matrix theory, whereas 
quantum analogs of classically integrable systems display the characteristics 
of the Poisson statistics 9 )~ 12 ). Some results in this direction for field theories 
have been obtained by Halasz and Verbaarschot: they studied the QCD 
lattice spectra for staggered fermions and its connection to random matrix 
theory 13 ^. 

In this paper we study quantum chaos in a field-theory schematic model. 
We analyze the energy fluctuation properties of the spatially homogeneous 
SU(2) Yang-Mills-Higgs (YMH) system (see Ref. 1-4). We show that these 
fluctuations give a clear quantum signature of the classical chaos-order tran- 
sition of the system. 

The Lagrangian density of the SU(2) YMH system 14 ) is given by 

L = l(D^(D^) - V(cf>) - \F; v F^ a , (1) 
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where 

(D^^d^-igA^, (2) 

F; v = d,A a v - d v Al + ge ahc A\Al , (3) 

with T b = a b /2, b = 1,2,3, generators of the SU(2) àlgebra, and where the 
potential of the scalar field (the Higgs field) is 

\/(0)=/i 2 |0| 2 + A|0| 4 . (4) 

We work in the (2+l)-dimensional Minkowski space (/i = 0, 1, 2) and choose 
spatially homogeneous Yang-Mills and the Higgs fields 

d t A* = d l( j) = 0, i = l,2 (5) 

i.e. we consider the system in the region in which space fluctuations of fields 
are negligible compared to their time fluctuations. 

In the gauge A$ = and using the real triplet representat ion for the Higgs 
field we obtain 

L = f + \{A 1 + A 2 )- g 2 [\A\A\ - 1{A X ■ À 2 f+ 

+ (Al + A\)^ - (À 1 ■ 0) 2 - (À 2 ■ 0) 2 ] - V(í) , (6) 

where $= 2 , 3 ), À 1 = (A\, Af, A\) and À 2 = (A\, A 2 2 , A\). 

When /i 2 > the potential V has a minimum at |0| = 0, but for /i 2 < 
the minimum is at 




which is the non zero Higgs vacuum. This vacuum is degenerate and after 
spontaneous symmetry breaking the physical vacuum can be chosen O = 
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(0, 0, v). If A\ — qi, A\ — q 2 and the other components of the Yang-Mills 
fields are zero, in the Higgs vacuum the Hamiltonian of the system reads 

H = \{p{ + PÍ) + 9 2 v\q\ + q\) + \g\lql , (7) 

where pi = q\ and p 2 = q 2 - Here w 2 = 2g 2 v 2 is the mass term of the 
Yang-Mills fields. This YMH Hamiltonian is a toy model for classical non- 
linear dynamics, with the attractive feature that the model emerges from 
particle physics. In the next sections we analyze first the classical chaos- 
order transition of the YMH system and then its connection to the quantal 
fructuations of the energy levels. 

2 Classical chaos-order transition 

A classical chaos-order transition for the YMH system has been observed pre- 
viously by different authors: Savvidy used the Chirikov'criterion 1 ), Kawabe 
and Ohta studied the Lyapunov exponents 3 -* and Salasnich analyzed the 
quantal overlapping resonances 4 ^ . In this paper we study the chaotic be- 
haviour of this YMH system by using the Gaussian curvature criterion of the 
potential energy 16 ) and the Poincarè Sections 17 ) . 

At low energy the motion near the minimum of the potential 

V(q 1 ,q 2 )=g 2 v 2 (q 2 1 +q 2 2 ) + ±g 2 q 2 q 2 2 , (8) 

where the Gaussian curvature is positive, is periòdic or quasiperiodic and is 
separated from the instability region by a line of zero curvature; if the energy 
is increased, the system will be for some initial conditions in a region of 
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negative curvature, where the motion is chaotic. According to this scenario, 
the energy E c of chaos-order transition is equal to the minimum value of the 
line of zero gaussian curvature K(q ± , q 2 ) on the potential-energy surface. For 
our potential the gaussian curvature vanishes at the points that satisfy the 
equation 

%% - <5s» 2 - < 29V + + ™ - ^ - ° ■ < 9 » 

It is easy to show that the minimal energy on the zero-curvature line is given 
by: 

E c = V mm (K = 0,q l )=Qg 2 v\ (10) 

and by inverting this equation we obtain v c = (E /Qg 2 ) 1 ^ 4 . We conclude that 
there is a order-chaos transition by increasing the energy E of the system 
and a chaos-order transition by increasing the value v of the Higgs field in 
the vacuum (see also Ref. 2). Thus, there is only one transition regulated 
by the unique parameter E/(g 2 v 4 ). 

It is important to point out that in general the curvature criterion guar- 
antees only a local instability 16 ^ and should therefore be combined with the 
Poincarè sections 17 -* (see Ref. 18). The classical equations of motion of the 
YMH system are 

9i = Pi , <?2 = P2 , Pi = -2g 2 v 2 q 1 - g 2 q x q\ , p 2 = -2g 2 v 2 q 2 - g 2 qfq 2 ■ 

(11) 

We use a fourth-order Runge-Kutta method 19 ) to compute the classical tra- 
jectòries. The conservation of energy restricts any trajectory of the four- 
dimensional phase space to a three-dimensional energy shell. At a particular 



energy the restriction q\ = defines a two-dimensional surface in the phase 
space, which is called Poincarè section. Each time a particular trajectory 
passes through the surface a point is plotted at the position of intersection 
(Q2,P2)- We employ a first-order interpolation process to reduce the inaccu- 



In Figure 1 we plot the Poincarè sections for different vàlues of the Higgs 
vacuum v but with the same energy E and interaction g. Chaotic regions 
on the surface of section are characterized by a set of randomly distributed 
points and regular regions by dotted or sòlid curves. The pictures show 
that the parameter v plays an important role: for large vàlues it makes 
the system regular. In fact, if we increase the harmònic part of the YMH 
potential the effect of the nonlinear term becomes less important. These 
numerical calculations confirm the analytical predictions of the curvature 
criterion: with E — 10 and g — 1 we get the critical value of the onset of 
chaos v c = (E/Qg 2 ) 1 l A ~ 1.14, in very good agreement with the Poincarè 



3 Quantum signature of the chaos— order tran- 



In quantum mechanics the generalized coordinates of the YMH system satisfy 
the usual commutation rules [qk,Pi] = i$ki-> with k, l = 1,2. Introducing the 
creation and destruction operators 



racies due to the use of a finite step length 17 ). 



sections. 



sition 




(12) 
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the quantum YMH Hamiltonian can be written 15 ) 

H = H + l -g 2 V , (13) 

where 

H = ujíyàlài + à\a 2 + 1) , (14) 
V = ^- 2 {à l + àtf{à 2 + àt)\ (15) 

with u 2 = 2g 2 v 2 and [àk, àf] = Ski, k,l — 1,2. 

The most used quantity to study the local fluctuations of the energy 
levels is the spectral statistics P(s). P(s) is the distribution of nearest- 
neighbour spacings = {E i+1 — Ei) of the unfolded levels É^. It is obtained 
by accumulating the number of spacings that lie within the bin (s, s + As) 
and then normalizing P(s) to unity 9 -* -12 ). 

For quantum systems whose classical analogs are integrable, P(s) is ex- 
pected to follow the Poisson limit, i.e. P(s) = exp (— s). On the other 
hand, quantal analogs of chaotic systems exhibit the spectral properties 
of GOE with P(s) = (7r/2)sexp (— |s 2 ), which is the so-called Wigner 
distribution 9 ^" 12 ). The distribution P(s) is the best spectral statistics to 
analyze shorter series of energy levels and the intermediate regions between 
order and chaos. 

Seligman, Verbaarschot and Zirnbauer 20 ) analyzed a class of two-dimensional 
anharmonic oscillators with polynomial perturbation by using the Brody 
distribution 21 ^ 

P(s,to) = a(uj + l)s w exp(-as w+1 ) , (16) 
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with 

a = (Ti^r 1 • (17) 

üü -+- 1 

This distribution interpolates between the Poisson distribution (uo = 0) of 
integrable systems and the Wigner distribution (cu — 1) of chaotic ones, and 
thus the parameter uj can be used as a simple quantitative measure of the 
degree of chaoticity. 

We compute the energy levels {Ei} with a numerical diagonalization of 
the truncated matrix of the quantum YMH Hamiltonian in the basis of the 
harmònic oscillators 22 ^. If \n\n 2 > is the basis of the occupation numbers of 
the two harmònic oscillators, the matrix elements are 

< nÍ7i2|#oKn 2 >= w(m + n 2 + l^nA'^ , (18) 

and 

< n x n 2 \V\n x n 2 >= -L·[^n 1 {n 1 - 1)^^-2+^ + + 2 K' 1 n 1+ 2+( 2n i+ 1 K' 1 n 1 } x 



x [^n 2 (n 2 - l)5 <n2 _ 2 + V(n 2 + l)(n 2 + 2)^„ 2+2 + (2n 2 + 1)^„J . (19) 

The symmetry of the potential enables us to split the Hamiltonian matrix 
into 4 sub-matrices reducing the computer storage required. These sub- 
matrices are related to the parity of the two occupation numbers n\ and 
n 2 : even-even, odd-odd, even-odd, odd-even. The numerical energy levels 
depend on the dimension of the truncated matrix: we compute the numerical 
levels in double precision increasing the matrix dimension until the first 100 
levels converge within 8 digits (matrix dimension 1156 x 1156) 22 ^ 23 ^. 

We use the first 100 energy levels of the 4 sub-matrices to calculate the 
P(s) distribution. In order to remove the secular variation of the level density 



as a function of the energy E, for each value of the coupling constant the cor- 
responding spectrum is mapped, by a numerical procedure described in Ref. 
24, into one which has a constant level density: {Ei} — > {Ei} (unfolding pro- 
cedure). We use the following Standard procedure to avoid mixing between 
states of different symmetry classes: 1) the diagonalization is performed for 
each sub-matrix (first 100 levels for each sub-matrix) ; 2) the unfolding is done 
for each sub-matrix; 3) the spacings are calculated for each sub-matrix; 4) the 
spacings of the 4 sub-matrices are accumulated to plot the P(s) distribution. 

In Figure 2 we plot the P(s) distribution for different vàlues of the pa- 
rameter v. The figure shows a Wigner-Poisson transition by increasing the 
value v of the Higgs field in the vacuum. By using the P(s) distribution and 
the Brody function it is possible to give a quantitative measure of the degree 
of quantal chaoticity of the system. Our numerical calculations show clearly 
the quantum chaos-order transition and its connection to the classical one. 

4 Conclusions 

The chaotic behaviour of an homogenous YMH system has been studied both 
in classical and quantum mechanics. The Gaussian curvature criterion and 
the Poincarè sections show that the chaotic behaviour is regulated by the 
unique parameter E/(g 2 v A ). The YMH system has a order-chaos transition 
by increasing the energy E and a chaos-order transition by increasing the 
value v of the Higgs field in the vacuum. 

The nearest-neighbour spacing distribution of the energy levels confirms 
with great accuracy the classical chaos-order transition of the YMH system. 
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In particular, the Brody function shows a Wigner-Poisson transition for the 
P(s) distribution in correspondence to the classical chaos-order transition. 

We observe that, as stressed previously, our YMH system is a toy model 
but it is very useful because it is possible to compare classical to quantum 
chaos. In the future will be important to study classical and quantum chaos 
in more realistic field theories. 
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Figure Captions 



Figure 1: The Poincarè sections of the model. From the top: v — 1, v — 1.1 
and v = 1.2. Energy E = 10 and interaction g = 1. 

Figure 2: P(s) distribution. From the top: v — 1 (u — 0.92), f = 1.1 
(üú = 0.34) and v = 1.2 (a; = 0.01), where u is the Brody parameter. First 
100 energy levels and interaction g — 1. The dotted, dashed and sòlid curves 
stand for Wigner, Poisson and Brody distributions, respectively. 
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